Dimerization induced by the RKKY interaction 
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We report the presence of spin dimerization in the ground state of the one dimensional Kondo 
lattice model at quarter filling. The emergence of this new phase of the Kondo lattice can be 
traced to the form of the RKKY interaction between the localized moments and provides the first 
example of dimerization induced indirectly by itinerant electrons. We propose this dimer ordering 
as the driving mechanism of the spin-Peierls phase observed in the quasi-one-dimensional organic 
compounds (Per)2M(mnt)2 (M=Pt, Pd). Moreover, this suggests that a richer phase diagram than 
the Doniach paradigm may be needed to accommodate the physics of heavy fermion materials. 
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Our understanding of the physics of heavy fermion ma- 
terials has been strongly influenced by the paradigm set 
by Doniach It proposes that the ground state of 
the Kondo lattice model exhibits, as a function of the 
Kondo coupling constant, antiferromagnetic and param- 
agnetic phases separated by a quantum critical point. A 
great deal of effort, both experimental and theoretical, 
has been devoted to the elucidation of the nature of this 
quantum phase transition (for an overview, see Q)- On 
the other hand, the variety of behavior observed in real 
materials suggests that the complete picture is far richer. 

We illustrate this richness by identifying a new phase 
of the one-dimensional Kondo lattice model Our nu- 
merical calculations with the density matrix renormal- 
ization group (DMRG) technique show that the ground 
state of the model has spin dimerization at quarter flll- 
ing (n = 1/2). We explain how this previously overlooked 
state can be understood in terms of the effective spin-spin 
RKKY interaction mediated by the conduction elec- 
trons. Interestingly, the quasi-one-dimensional organic 
compounds (Per)2M(mnt)2 (M=Pt, Pd) are actual re- 
alizations of weakly coupled quarter-filled S = ^ Kondo 
chains (EBi00|- We will show, based on our results, 
that the itinerant-electron-induced dimer order we have 
found is a viable candidate for the driving mechanism 
behind the hitherto unexplained dimerization transition 
observed in these systems. 
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where Cja annihilates a conduction electron in site j with 
spin projection a, Sj is a localized spin-i operator and 
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Cape, g is the conduction electron spin 
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density operator. J > is the Kondo coupling constant 
between the conduction electrons and the local moments. 
We have set the hopping amplitude and the lattice spac- 
ing to unity to fix the energy and length scales. We 



treated the model with the DMRG technique with 
open boundary conditions. We used the finite-size algo- 
rithm for sizes up to L = 120 keeping up to m = 800 
states per block. The discarded weight was typically 
about 10^^ — 10^^ in the final sweep. 
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Figure 1: (a) The dimer order parameter D{j) vs lattice site 
for J = 0.5, L = 120 and density n = 0.5. Only half the chain 
is shown, (b) D{L/2 - 1) vs 1/L for J = 0.5 and n = 0.5. 



Dimerization of the localized spins may be detected 
Sj (1) through the order parameter D [j) = (Sj • Sj+i). In a 



uniform system D{j) is j-independent, whereas it shows 
oscillations of period 2 in the presence of dimerization. 
In Fig. n](a) we show D{j) at quarter filling {n = 1/2) 
for J = 0.5 and L = 120. There are strong oscillations 
with amplitude « 0.21 around an average zero value. We 
have checked that this is not an artifact of a finite system 
or open boundaries. Indeed, in Fig. n](b), we show the 
value of the order parameter at the center of the chain 
D{j = L/2—1) as a function of 1/L. There is very weak 
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size dependence indicating that the value aX L = 120 
is already very close to the thermodynamic limit. This 
robust result is unambiguous evidence for a dimerized 
ground state at quarter-filling and it is surprising that 
it has gone unnoticed in such a well studied model 
As far as we know, it is the first example of dimerization 
induced by an indirect spin-spin interaction mediated by 
itinerant electrons. 
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Figure 2: (a) Modulus of the smoothed Fourier transform 
|-D(g)| vs momentum for J = 0.5, L = 120 and density n = 
0.5. (b) \D{q = 7r)| as a function of the Kondo coupling J for 
L = 120 and n = 0.5. 

It is also instructive to study the Fourier transform 
of D{j). In order to minimize boundary effects we 
used a smoothed Fourier transform as suggested in Refs. 
[H m. In Fig. 12 (a) we show the smoothed Fourier 
transform of D[j) for J = 0.5, n = 0.5 and L = 120. 
The only feature is a noticeable peak at g = tt refiecting 
the alternating sign of the order parameter D{j). The 
critical value of J that separates the paramagnetic phase 
from the ferromagnetic one is Jc ~ 1.7 for n — 0.5 
We found that the dimerization occurs only inside the 
paramagnetic phase. In Fig. |21 (b) we show the inten- 
sity of the peak of D{q) at g = tt as a function of the 
Kondo coupling J iov L = 120 and n — 0.5. As can be 
seen, the dimerization decreases smoothly as we increase 
J and disappears in the ferromagnetic phase (within the 
accuracy of the DMRG). 

A very small magnetic field is enough to kill the dimer- 
ization. We show in Fig.^lthe magnetization density and 
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Figure 3: Magnetization density M and modulus of the 
smoothed Fourier transform \D{q = 7r)| as a function of ap- 
plied magnetic field h, for L = 40, J = 0.5 and n — 0.5. Inset: 
Magnetization density up to higher fields. 



the smoothed Fourier transform of the dimer order pa- 
rameter as a function of magnetic field h for J = 0.5 
and n = 0.5 (making g^iB — I for both localized spins 
and conduction electrons). Since finite-size effects are 
negligible, the results here are for L = 40. There is a 
characteristic field h* « 0.05, which separates a high 
susceptibility region at h < h* from a low susceptibil- 
ity one at h > h* (see inset of Fig. EJ. Similar results 
were observed for other densities and Kondo couplings. 
As h is increased up to h*, the magnetization density 
grows linearly up to M « 0.5. This corresponds to the 
full polarization of the localized spins. A further increase 
of h past h* acts to slowly polarize the conduction elec- 
trons (eventually reaching full saturation at M = 0.75), 
which are now effectively decoupled from the local mo- 
ments. This is confirmed by the value of the slope for 
h h* , which is very close to the spin susceptibility of 
the non-interacting one-dimensional electron gas. Thus, 
the field h* is the characteristic energy scale of the in- 
teraction between the spin lattice and the electron sea. 
The numerical values of h* and the ratio of high-field 
and low-field susceptibilities (« 0.008) both suggest the 
presence of an exponentially small scale, as seen in other 
studies 01 ■ This scale also governs the suppression of the 
dimerization as can be seen in Fig.^l Since at h> h* the 
conduction electrons effectively decouple from the local- 
ized spins, this is direct evidence that the dimerization is 
induced by the interaction with the itinerant electrons. 

Dimerization at quarter-filling is perhaps not too sur- 
prising if we consider the RKKY interaction between 
localized spins a distance j apart. This is appropriate 
if J <g; 1. Using the usual RKKY formula for a one- 
dimensional tight-binding lattice at quarter fiUing we can 
get the first ten couplings Jrkky (j) as shown in Fig. 01 
Focusing initially on the first two neighbors, we see that 
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Figure 4: RKKY coupling constant (in arbitrary units) as a 
function of the distance between sites. 



Figure 5: Modulus of the smoothed Fourier transform D{q) 
away from quarter filling {L = 120 and J = 0.35). 



J RKKY (1) < and J RKKY (2) > 0. It is believed that a 
spin chain with first and second neighbor interactions 
Ji and J2 is spontaneously dimerized if J2 > and 
-4J2 < Ji < and this is satisfied by the RKKY 
couplings. In fact, even if Ji > dimerization would 
ensue provided that Ji < 4.15 J2 0,0- 

Another curious aspect of Fig.^is the sign alternation 
of D{j) with an almost constant amplitude ~ 0.21. A 
classical localized spin configuration of the type TTiiTT 
• • • would give rise to such an alternation with ampli- 
tude 0.25. Interestingly, the further neighbor RKKY 
couplings of Fig. 31 favor this classical state (note that 
J RKKY (j) ~ for odd j > 3). Of course, we do not 
expect this long-range classical order to survive in the 
singlet ground state, but the fact that |-D(j)| ^ 0.25 sug- 
gests that it may have a very long correlation length. 

We have also calculated dimer correlations away from 
n = 0.5. We observed that D(j) no longer oscillates sym- 
metrically around zero. The average actually goes from 
positive (ferromagnetic) for n < 0.5 to negative (anti- 
ferromagnetic) for n > 0.5. This is consistent with the 
general trend of the RKKY interaction which is predomi- 
nantly ferromagnetic at low fillings. However, in contrast 
to the quarter-filled case, away from n = 0.5 D{j) has 
more than one characteristic wave vector. This is shown 
in Fig. for the densities n = 0.4 and n = 0.6 (L = 120 
and J = 0.35). They have Akp — 27rn (mod 27r) as the 
dominant peak. There are also sub-dominant peaks at 
2i^smaU _ 2k p^^'^ = 7r(l — n), which are the Fermi 

sea sizes without and with the locaHzed moments in- 
cluded in the count, respectively 0,0,0]. However, no 
clear trend can be discerned. Quantum fiuctuations will 
destroy incommensurate order in one dimension. How- 
ever, the system may lock onto particular commensurate 
structures within finite intervals of n. Such long com- 
mensurate periods have been discussed in related Kondo 
lattice models 0| used to describe the pnictide CeSb jlDl 
and in the context of doped Mott insulators jiJi • 

The organic compounds (Per)2M(mnt)2, where M=Pt 



or Pd, "Per" is perylene and "mnt" is maleoni- 
triledithiolate, are realizations of quarter-filled quasi-one- 
dimensional Kondo lattices 0,0, Q. The locaHzed spins 
reside in the M(mnt)2 units and the conduction electrons 
are provided by the perylene groups. These compounds 
show metal insulator transitions at Tmit (Pt) — IK 
and Tmit (Pd) = 28 K, accompanied by one-dimensional 
lattice instabilities which indicate the formation of a 
Peierls-type dimerization. Initial studies jl| had difficul- 
ties reconciling the behavior of these systems with con- 
ventional electron-phonon (Peierls) or spin-phonon (spin- 
Peierls) instabilities and the exact nature of these transi- 
tions has remained a mystery. On the other hand, there 
were speculations on the role played by the RKKY in- 
teraction 01 . Several uncertainties remained, however, 
since spin-density wave (SDW) ordering was not observed 
and the Fermi wave vector of the conduction electrons 
{2kp = 1:12) does not match the period of the dimeriza- 
tion. 
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Figure 6: Charge gap vs 1/L for various fillings. The notation 
is the same as in ref. The accuracy is given by the size 

of the symbols. 

Our results suggest a possible way out of this quandary. 
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Indeed, the one-dimensional nature of the system sup- 
presses the appearance of a conventional SDW and favors 
a dimerized state instead. Furthermore, the Fermi wave 
vector, as explained in connection with Fig. ^ is actu- 
ally crucial for providing the correct spatial dependence 
of the RKKY interactions that lead to the dimer order. 
An estimate of the finite transition temperature (absent 
in the strictly one-dimensional model) could be obtained 
from the divergence of the dimer susceptibility as T 0, 
along the lines explained in Refs. (22, 1231 ■ Our results 
on the ground state, however, suggest the presence of 
a (probably exponentially) small energy scale of a size 
compatible with the transition temperatures (Fig.EJ. It 
is therefore tempting to use this scale as an estimate of 
Tmit- Moreover, the direct measurement of the mag- 
netic field needed to destroy the dimer order at T = 
showed it to be of the order of ksTMiT / 0, lending 
further support to this identification. 

Finally, in Fig. we plot the finite size charge gap for 
several densities (in the sector of total spin zero) as a 
function of It is clear that the quarter-filled case 

has a charge gap, unlike other fillings (see Ref. |0| for 
a similar case). A mechanism for the opening of this 
gap can be obtained by bosonization. By analogy with 
the RKKY interaction, if we integrate out the localized 
spins an effective interaction between conduction elec- 
tron spins at sites i and j is generated, which is pro- 
portional to (Si • Sj) . The strong dimer correlations of 
Fig. n lead to a term (x ^ij+i (—!)■'. This gives rise 
to a staggered nearest neighbor interaction between con- 
duction electron spins (—1)"' Sj ■ Sj+i. Bosonizing this at 
quarter-filling leads to a cos (\/8Kp(j)p) term, where we 
follow the notation of Ref. ,24]- If Kp < 1 this term is 
relevant and opens a gap in the charge sector in close 
analog y to the umklapp term in the half-filled Hubbard 
model \M. This is a very mild condition and is likely to 
be fulfilled. Indeed, numerical estimates of Kp at other 
fillings indicate that Kp < 1 Thus, we propose 

this RKKY-induced dimerization as the driving mecha- 
nism behind the metal-insulator transitions observed in 
(Per)2M(mnt)2. Note that a similar term is generated 
in the spin sector cos [^/SK^cfia) ■ However, Ka — 1 due 
to SU(2) symmetry and this is a marginal term. Its rel- 
evance and the accompanying spin gap depends on the 
sign of the coupling constant. We have found no evidence 
of a spin gap at this filling. 

In conclusion, we have found that the Kondo lattice 
model in one dimension has a novel type of dimeriza- 
tion at quarter conduction electron filling. It can be un- 
derstood from the structure of the RKKY interaction 
between localized moments. Furthermore, it provides 
a simple mechanism for the metal-insulator transition 



of some quasi-one-dimensional organic compounds. Al- 
though confined to one dimension, we believe our results 
may have implications for the phase diagram of three- 
dimensional heavy fermion materials. 
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